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Magnetic skyrmions are localized, topologically protected spin-structures that have been proposed for stor¬ 
ing or processing information due to their intriguing dynamical and transport properties. Important in terms 
of applications is the recent discovery of interface stabilized skyrmions as evidenced in ultra-thin transition- 
metal films. However, so far only skyrmions at interfaces with a single atomic layer of a magnetic material 
were reported, which greatly limits their potential for application in devices. Here, we predict the emergence 
of skyrmions in [4d/Fe2/5d]n multilayers, i.e. structures composed of Fe biatomic layers sandwiched between 
4d- and 5d-transition-metal layers. In these composite structures, the exchange and the Dzyaloshinskii-Moriya 
interactions which control skyrmion formation can be tuned separately by the two interfaces. This allows engi¬ 
neering skyrmions as shown based on density functional theory and spin dynamics simulations. 


Topologically protected spin-structures, in particular mag¬ 
netic skyrmions CH3, have recently received much atten¬ 
tion due to their promising applications in spintronics E- 
m. Among the different mechanisms that can induce a 
topologically non-trivial magnetic structure, the relativistic 
Dzyaloshinskii-Moriya interaction (DMI) (Tl O arising in 
non-centrosymmetric systems is maybe the most prominent 
one 13. Magnetic skyrmions have been first revealed in a 
few non-centrosymmetric bulk magnets, such as the chiral 
MnSi |[T0l[TTl or FeCoSi |[T2l. In two-dimensional (2D) sys¬ 
tems the DMI arises always due to the broken inversion sym¬ 
metry at the interface between a magnetic film and a sub¬ 
strate CMa. E.g. an Fe monolayer on an Ir(lll) substrate 
exhibits strong DMI and on this system a nanoscale skyrmion 
lattice was first observed by scanning-tunneling microscopy 
(STM) ESI. This skyrmion lattice is stabilized by the short- 
range four-spin interaction and slight modifications of the film 
composition, e.g. adding an atomic Pd overlayer, can push 
this system into a regime where individual skyrmions can ex¬ 
ist which can be manipulated by spin-polarized currents in 
STM 03. 

The manipulation of ultrathin films such as Fe/Ir(lll) is 
a versatile route to tailor properties such as the size or mag¬ 
netic field needed to induce a magnetic skyrmion 
For practical applications, however, there are several difficul¬ 
ties that have to be overcome: (i) Manipulation of skyrmions 
by electric currents El EMI and measurement of the topo¬ 
logical Hall effect ll25HZ7l requires that a substantial fraction 
of the current runs through the magnetic structure. In a sin¬ 
gle Fe layer on a metallic substrate this is difficult to realize, 
(ii) Creating large area single Fe layers on a substrate requires 
specialized preparation methods that hamper large-scale pro¬ 
duction of these structures, (iii) Temperature stability and op¬ 
timization of the size of these magnetic structures requires a 
tuning of different parameters such as exchange-interactions, 
the DMI, and the magnetic anisotropy. 

Here, we propose multilayers composed of Fe bilayers epi¬ 
taxially sandwiched between Ad- and 5 ti-transition-metal lay¬ 
ers as promising systems towards solving these issues. The 


repetition of layered structures is a powerful method in spin¬ 
tronics to engineer materials properties (281 [29l: Besides in¬ 
creasing the magnetic material in [Ad/Fe/5d]n multilayers, 
such structures have the key advantage that the essential mag¬ 
netic interactions, the exchange and the DMI, can be con¬ 
trolled by two separate interfaces. Thereby, the requirement 
of pseudomorphic growth is limited to the interface which in¬ 
duces the DMI while intermixing at the other interface is not 
crucial. Based on density functional theory (DFT), we show 
a route to engineer the skyrmion systems by varying the in¬ 
terface composition. We focus on an Fe bilayer sandwiched 
between Ir and Rh/Pd layers, however, other combinations of 
Ad and bd transition-metals are also feasible. Using spin dy¬ 
namics simulations we find the emergence of skyrmions in the 
Fe bilayers and study their properties. 

In order to obtain the magnetic interactions in transition- 
metal (TM) multilayers we have performed DFT calculations 
for various noncollinear spin structures including spin-orbit 
coupling (301 EH- By mapping the total energies of the differ¬ 
ent magnetic configurations to an extended Heisenberg model 
we can extract parameters for the interactions. We consider 
homogeneous flat spin spirals with a period A in which the 
magnetization rotates along a given crystallographic direc¬ 
tion that lies in the plane of the interface by an angle of 
(f = {27r/X)d from atom to atom where d is the spacing be¬ 
tween adjacent atoms (see methods). This approach has been 
successfully used in the past to study complex spin textures 
in a wide range of systems such as biatomic chains at sur¬ 
faces (32l, ultra-thin film skyrmion systems (161 [181 and mag¬ 
netic mono- |[T5l[33l[3I and bilayers (35l. 

A characteristic quantity of the exchange interaction in the 
considered TM multilayers is the effective nearest-neighbor 
exchange constant Jeff. It is a measure of how quickly the en¬ 
ergy rises for long-wave length spin spirals with respect to the 
ferromagnetic state. Jeff was obtained by fitting the dispersion 
curve ^(q) of spin spirals characterized by the wave vector 
q in the vicinity of the ferromagnetic (FM) state (|q| = 0) 
with a parabola Jeff^|gp where a is the nearest-neighbor 
distance within the layers (for energy dispersions see supple- 
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mentary material S2). Therefore, a negative value of Jeff in¬ 
dicates a spin spiral ground state driven by exchange interac¬ 
tion. As a second characteristic quantity, we include the en¬ 
ergy difference between a spin spiral state with a wave length 
X = 27r/\q\ = 2.7 nm and the ferromagnetic state, denoted by 
^ss- 

Jeff and J^ss are displayed for various systems in Fig. 
We start with the monolayer films on substrates which can 
host single skyrmions as shown in experiment ifTTl and ex¬ 
plained by theory d). For both fee and hep stacking of the 
Pd overlayer on Fe/Ir(lll), we find a small value of the ef¬ 
fective exchange constant, i.e. Jeff = —2.3 and +4.4 meV, 
respectively. The DMI which is of similar strength in both 
systems lIT^ can stabilize skyrmions in these films in an ap¬ 
plied magnetic field. For two layers of Pd on Fe/Ir(l 11), Jeff is 
negative and slightly stronger than for Pd(fcc)/Fe/Ir(lll) and 
a skyrmion phase also occurs in a magnetic field as found in 
our simulations. In these two systems, the spin spiral ground 
state is already enforced by the exchange interaction and be¬ 
comes even more favorable due to the DMI. We can also con¬ 
clude the spin spiral ground state from the negative value of 
Ess which closely follows the trend of Jeff. 
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FIG. 1: Effective exchange interaction Jeff (full circles) in transition- 
metal multilayers (black/blue fonts) and ultra-thin films on surfaces 
(gray fonts). Insets show sketches of the multilayer and thin film 
geometry. The energy difference between a spin spiral with a wave 
length of 2.7 nm with respect to the ferromagnetic state, J^ss is given 
by empty squares. The dashed line is a guide for the eye. 


Now we turn to the multilayers which are built by an epi¬ 
taxial fee stacking of (111) layers with a geometry as in an 
Ir bulk crystal but allowing for structural relaxation in the 
direction normal to the layers (see methods). We first con¬ 
sider a multilayer which is a repetition of a sandwich con¬ 


sisting of one atomic Fe layer between two Pd and two Ir 
layers denoted as 2Pd/Fe/2Ir. As shown in Fig. the effec¬ 
tive exchange interaction is very similar to that of the ultra- 
thin film system Pd(fcc)/Fe/Ir(lll) which demonstrates that 
Jeff is mainly controlled by the interfaces. A first naive so¬ 
lution to increase the amount of magnetic material, i.e. Fe 
layers in the sandwich, would be to repeat the Pd/Fe inter¬ 
face sequence as in Pd/Fe/Pd/Fe/2Ir. However, Jeff becomes 
much larger (?^ 15 meV) indicating strong ferromagnetic ex¬ 
change interaction which prevents skyrmion formation (see 
supplement S4). If on the other hand, we consider the mul¬ 
tilayer with a basic unit consisting of 2 layers of Pd, Fe, and 
Ir, i.e. 2Pd/2Fe/2Ir, the effective exchange interaction is of the 
same order as in the original single Fe-layer system. 

We can understand the role of Pd by replacing it by Rh 
which has one electron less in the 4J shell. For a multi¬ 
layer consisting of 2Rh/2Fe/2Ir we obtain a value of Jeff ~ 
— 16 meV, i.e. negative and of large absolute value. In this 
system, we find an exchange driven spin spiral ground state 
and no skyrmion formation in a reasonable external magnetic 
field. By choosing an intermediate interface composition, 
e.g. Rh/Rho.75Pdo.25/2Fe/2Ir, we can increase the Ad band fill¬ 
ing quasi continuously and thereby increase the value of the 
effective exchange coupling Jeff. For a large range of inter¬ 
face alloy compositions Rh^^Pdi-a^, we find multilayers with 
a small Jeff which can host skyrmionic ground states as we 
show below using spin dynamics simulations. 

After having revealed the feasibility to tune the magnetic 
interactions in transition-metal multilayers, we demonstrate 
the formation of skyrmions in such systems. In order to ex¬ 
plore the magnetic phase space in an external magnetic field, 
we employ the spin Hamiltonian 

^ = - E ■ M,) - ^D,, ■ {Mi X M,) 

ij ij 

- Y.K{Mtf ( 1 ) 

i i 

which describes the magnetic interactions between the mag¬ 
netic moments of atoms at sites R^. The parameters 
for the exchange interaction (J^j), the DMI (D^^) as well as 
an uniaxial magnetocrystalline anisotropy (AT) were obtained 
from DFT. 

In a magnetic bilayer it is convenient to split the exchange 
interaction into an intralayer contribution, j|, which corre¬ 
sponds to the exchange coupling to an atom of the 6-th near¬ 
est neighbor within an Fe layer, and the interlayer part J^, 
which corresponds to the coupling between the two Fe layers 
(see supplementary material S2 for details). For the multi¬ 
layer Rh/Pd/2Fe/2Ir, on which we focus below as a generic 
example, the interlayer coupling between nearest neighbors 
is very strong J^- = 24.7 meV and favors ferromagnetic 
alignment of magnetic moments of the two Fe layers. How¬ 
ever, the nearest-neighbor exchange coupling within each Fe 
layer, j| = 3.9 meV, is much smaller and on the order 
of antiferromagnetic exchange with third nearest neighbors. 
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jI = —3.0 meV. In total, the competition between intra- and 
interlayer exchange leads to the small effective exchange con¬ 
stant Jeff = 1 meV (see Fig. [T]), which indicates a very flat 
dispersion curve close to the FM state (cf. supplement S3). 

The DMI is of primary importance to stabilize skyrmions 
and has to be carefully evaluated for the bilayer case. Al¬ 
though there are contributions from the two interfaces, we And 
from our DFT calculations that the DMI is controlled by the 
Fe interface to the heavy 5J material (see supplement S3). 
Within our spin Hamiltonian, we can describe it by an effec¬ 
tive nearest-neighbor D-vector acting on the Fe moments at 
the Fe/Ir interface only. We obtain a value of D = 1.3 meV 
for Rh/Pd/2Fe/2Ir with very little variation upon changing 
the composition of the Rh^^Pdi-a^ interface (on the order of 
0.1 meV). This value is very similar to those obtained for the 
ultra-thin film system Pd/Fe/Ir(lll), i.e. D = 1.0 meV and 
1.35 meV for fee and hep stacking of Pd, respectively (m. 
The strength of the DMI is sufficient to create a spin spiral 
ground state for Rh/Pd/2Fe/2Ir with a right hand rotation sense 
and a period of A = 2.25 nm (cf. supplement S3). The ob¬ 
tained uniaxial anisotropy, K, is 0.6 meV per Fe atom and 
favors an out-of-plane magnetization. 



B(T) 

FIG. 2: Low temperature magnetic phase diagram of the 
Rh/Pd/2Fe/2Ir multilayer shown in the inset. The total energy of 
the ferromagnetic (FM) state (green line) and of the skyrmion lattice 
(SkX, red line) is shown with respect to the spin spiral (SS) state 
(dashed black line). The regimes of the SS, SkX, and FM phase are 
indicated by blue, red, and green color, respectively. 

We show the formation of skyrmions in multilayers 
by studying our spin Hamiltonian parametrized from first- 
principles using spin dynamics simulations (see methods). In 
Fig. 1^ the obtained low temperature phase diagram is shown 
as a function of an external magnetic field applied perpen¬ 
dicular to the layers. At 5 = 0 T the lowest energy is ob¬ 


tained for a spin spiral with a period of A = 2.25 nm. As 
the magnetic field increases, the spin spiral is destabilized 
and at B = 2.95 T a transition to a skyrmion lattice oc¬ 
curs. Fields of this strength have been used e.g. in ReffTTlto 
induce skyrmions in monolayers experimentally and further 
fine-tuning of the materials in the multilayer can be used to 
get to even lower critical fields. At B = 5.1 T, the magnetic 
field becomes large enough to stabilize the FM state. 

Since the DMI does not vary much upon changing the com¬ 
position X of the Rh/Rha^Pdi_a^/2Fe/2Ir multilayer and the in¬ 
fluence of the magnetocrystalline anisotropy energy is rather 
moderate, the magnetic field strengths at the phase transitions 
depend mainly on the effective exchange interaction Jeff. Al¬ 
though this interaction is influenced by factors such as alloy¬ 
ing at the interfaces, it can be tuned as shown in Fig. [fallow¬ 
ing a control over the phase diagram. (Limiting scenarios with 
respect to variation of the exchange interaction and magne¬ 
tocrystalline anisotropy are discussed in the supplement S4.) 
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FIG. 3: Magnetization profile of an isolated skyrmion in the multi¬ 
layer Rh/Pd/2Fe/2Ir averaged at a given radial distance r from the 
center. The polar angle 0 of the magnetization vector is shown as a 
function of r. Red dots correspond to the Fe atoms at the Fe/Ir in¬ 
terface and the black dots correspond to the Fe atoms on the Pd side. 
The inset shows a top view of the magnetization of the skyrmion with 
red arrows pointing up and green/blue arrows pointing down. 

In Fig. If we present the profile of an isolated skyrmion in 
the bilayer at B = 3.5 T, i.e. the polar angle 0 of the local 
magnetization as a function of the distance from the center. 
The curves for both Fe layers show that at the center of the 
skyrmion the magnetic moments are antiparallel to the sur¬ 
rounding FM background. We obtain a skyrmion diameter 
of approximately 2.2 nm consistent with the spin spiral min¬ 
imum. The similarity between the two curves indicates that 
the magnetic moments in the two layers are strongly coupled 
and that the skyrmion exists in both layers simultaneously. 

In order to understand the stabilization mechanism of a 
skyrmion in the Fe bilayer, we study the contributions of the 
different interactions to the energy density with respect to the 
ferromagnetic background in each Fe layer for the isolated 
skyrmion at B = 3.5 T. The energy density due to exchange 
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FIG. 4: Energy density contributions for an isolated skyrmion in an 
Rh/Pd/2Fe/2Ir multilayer at a magnetic field of B = 3.5 T. (a,b) 
show the exchange energy in the Fe layer at the Ir interface (Fe@Ir) 
and at the Pd interface (Fe@Pd), respectively. (c,d) show the total en¬ 
ergy (exchange, DMI, Zeeman and anisotropy) for the two Fe layers 
and (e) for the complete Fe bilayer. The color scale bar of the energy 
is given to the right of the panels indicating stabilizing {E < 0) and 
destabilizing {E > 0) contributions. For (e) the scale bar is the same 
as for (c) and (d). The insets show an average of the energy density 
over the 6^^ neighbor shell, (f) total energy density averaged over 
the neighbor shell for (e). The inset shows the integrated total 
energy of the Fe bilayer (red), Fe@Pd (green) and Fe@Ir (black). 


(Figs.|^ ,b) is very similar for the two Fe layers apart from a 
rotation by 60° due to the different position of the first near¬ 
est neighbor atoms in the adjacent layer which are coupled 
by giving the main contribution to the exchange. The 
exchange energy density is minimal in the skyrmion center 
where a moderate rotation of the magnetization (cf. Fig. 
corresponding to a long wavelength spin spiral is energetically 
favored. As we move away from the center of the skyrmion, 
the canting between nearest-neighbor spins and the exchange 
energy increases, reaches a maximum at a distance of 1 nm 
from the core and decreases again when reaching the FM 
background. 

The total energy density reveals the difference between the 
two Fe layers. For the Fe layer at the Pd interface (Fig.|^) the 
DM interaction vanishes, and the total energy density is only 
increased by the Zeeman term and the anisotropy contribution. 
However, the Fe layer at the Ir interface (Fig.|^) is subject to 
a significant DMI and its energy contribution outweighs the 
anisotropy and the Zeeman term leading to a stabilizing (neg¬ 


ative) energy density within the entire skyrmion. The total 
energy density of the Fe bilayer, i.e. the sum of the two lay¬ 
ers, shown in Figs. & and f, exhibits a small local minimum 
at the center of the skyrmion due to the compensation of the 
DMI from the Fe/Ir interface and the Zeeman and anisotropy 
terms. As we move away from the center, the total energy 
density first increases due to the strong exchange and Zeeman 
contributions and the rise of the anisotropy term and decreases 
beyond approximately 1 nm due to the DM term (cf. Fig.|^. 

The energy density integrated over a circle centered at the 
skyrmion is shown as a function of the radius r for the Fe bi¬ 
layer and the two Fe layers in the inset of Fig. |^. For the Fe 
layer at the Pd interface (green curve), the integrated energy 
is always positive as expected from the energy profile (inset 
of Fig. 1^) and saturates at approximately r = 1.25 nm. For 
the Fe layer at the Ir interface (black curve), the energy is 
negative indicating a skyrmion stabilization due to the DMI. 
The competition of these two contributions leads to the pro¬ 
file for the Fe bilayer (red curve) which displays a positive 
maximum at small radius but converges to a negative value 
of about —32 meV indicating that the skyrmion is energeti¬ 
cally favorable due to the DM contribution at the rim of the 
skyrmion. 

We have focused on the system composed of an Fe bilayer 
sandwiched between Rh/Pd and Ir layers. For Fe monolayers 
on M- and 5(i-transition metal surfaces we have demonstrated 
before that the exchange interaction can be tuned by the sub¬ 
strate band filling from ferro- to antiferromagnetic (SbHSHl. 
Therefore, other TM multilayers, e.g. Rua^Pdi_a^/2Fe/Pt allow 
a similar tuning of Jeff by varying the composition x. The 
strength of the DMI, on the other hand, will be determined by 
the interface with the 5 J-TM, e.g. Ir or Pt |[39l|40l. Transition- 
metal multilayers consisting of a repetition of a few atomic 
layers as proposed here have recently been realized experi¬ 
mentally to stabilize chiral domain walls 1^ . 

In conclusion, we have shown that magnetic skyrmions can 
emerge in transition-metal multilayers consisting of Fe bilay¬ 
ers sandwiched between 4J- and 5J-transition metals. Such 
systems provide a rich field for the formation of skyrmions 
and allow to go beyond the monolayer systems discussed so 
far EMSl to a thickness range where the unique transport 
properties of these systems can be studied and exploited in 
devices. Recent success in tailoring the DMI at the interface 
to the 5J material 1^ can be combined with our concept to 
tune the exchange interaction by the choice of the M com¬ 
ponent of the multilayers, thus opening a route to engineer 
skyrmion properties for spintronic applications. 


Methods 

We have explored the multilayer structures from first- 
principles based on the full-potential linearized augmented 
plane wave method as implemented in the FLEUR code IdTIl . 
Within this approach we can calculate the total energy of non- 
collinear magnetic structures such as spin spirals |[30l includ- 
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ing the DMI in first order perturbation theory with respect to 
the spin-orbit coupling ED. We have used a two-dimensional 
hexagonal p(l x 1) unit cell within each layer and the in plane- 
lattice parameter of the Ir(l 11) surface as obtained from DFT 
in Ref. HH). The distances between the different layers were 
relaxed using the mixed functional suggested in Ref. 14^ . 
which is ideally suited for interfaces of 3d- and 5d-transition 
metals. The magnetic properties were obtained within the lo¬ 
cal density approximation ||4^ . In order to study alloying 
within the RhPd layers, we have used the virtual crystal ap¬ 
proximation (VGA) 1441 which has been used successfully 
for various types of systems from solid solutions ma to mag¬ 
netic ultrathin films 1^ . Further computational details can 
be found in the supplementary material (SI). 

In order to calculate the different energy differences be¬ 
tween the FM, SkX and SS phases, we have relaxed a spin 
lattice of the bilayer with (lOOx 100) spins in each layer ac¬ 
cording to the Landau-Lifshitz equation of spin dynamics: 


M, ^ dHi 
dt dMi 


xMi- p 



X Mi 


( 2 ) 


where P is the damping term and Hi is the Hamiltonian of 
spin i and is expressed by the Hamiltonian given in Eq. (1) 
if we omit the sum over i. In order to achieve fast relax¬ 
ation, we have used (3 = 0.95. We have carried out sim¬ 
ulations on a micro-second time scale with time steps rang¬ 
ing from 0.1 ps to 10 ps . The equation of motion was in¬ 
tegrated with a simple Euler, a Heun, and two semi-implicit 
integrators Bbl leading to the same results. To ensure that we 
reached the ground state, we have simultaneously run Monte 
Carlo simulations with a Metropolis algorithm on these con¬ 
figurations at T = 0 K. Depending on the system, both tech¬ 
niques gives similar results with a precision on the order of a 
few yueV/spin, sufficient to discriminate between the different 
phases in our cases. 
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